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Abstract 
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1 Introduction 

The antiferromagnetic spin-^ XXZ chain with Hamiltonian 

7=1 

is a rare exception, even among the models solvable by Bethe ansatz and connected 
to the Yang-Baxter equation, in that much is known exactly about its correlation func- 
tions. All the remarkable progress achieved in the past decade originated in the deriva- 
tion [13, 14] of a so-called multiple integral representation for the density matrix of the 
model. The next key works were [21], where the results of [13, 14] were rederived by 
means of the algebraic Bethe ansatz and generalized to include a non-zero longitudinal 
magnetic field, then the article [5], where some of the multiple integrals were actually 
calculated for the first time, and [6], where an ansatz for the explicit form of a special 
density matrix element was formulated (for the inhomogeneous and isotropic model). 

The ideas developed in these works evolved further in various recent articles. The 
ansatz of [6] for the inhomogeneous generalization of the correlation functions was 
related to the g-Knizhnik-Zamolodchikov equation, it was generalized to the general 
density matrix element and finally proved (see [3,4] and the literature cited there). The 
technique of [5] was developed further, and a number of concrete examples of short 
range correlation functions of the XXZ chain was worked out, e.g. [15, 16]. Within the 
Bethe ansatz based approach integral representations for correlation functions other 
than the density matrix elements were developed [18], the important problem of the 
summation of density matrix elements into two-point functions was attacked success- 
fully [20], and multiple integral representations for time-dependent correlation func- 
tions were suggested [19]. 

Based on a combination of the Bethe ansatz approach with quantum transfer matrix 
techniques [1 1, 22, 23, 26] the authors in collaboration with A. Kliimper derived multi- 
ple integral representations for finite temperature correlation functions [9] and density 
matrix elements [8, 10]. 

This article is part of a programme of exploring the properties of finite temperature 
correlation functions from their multiple integral representations. So far the multiple 
integrals have proved to be useful for deriving high temperature expansions for correla- 
tion functions [28] and have also been evaluated numerically in the simplest non-trivial 
situation of three lattice sites [7]. Here we shall show how the multiple integrals re- 
duce to single integrals in the XX limit and in the Ising limit. As far as the XX limit 
is concerned our discussion in many respects parallels that for the zero temperature 
case [17]. 

The paper is organized as follows. In the next section we review the multiple 
integral formula for the one-parameter generating function [12] of the S z -S z correlation 
functions of the XXZ infinite chain, that was obtained in [9] and generalizes a formula 
in [18] to finite temperature. Then we treat the XX limit and the Ising limit in two 
separate sections. We close with a short summary. 
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2 Integral representation of the generating function 

Starting point of our considerations in the following two sections will be the multiple 
integral representation for a one-parameter generating function of the S z -S z correlation 
functions derived in [9] . This one-parameter generating function is the thermal average 

m 

4>(cp|m) = <exp{(p £ ejl}) Th , (2) 

7=1 

where e\ = (R ?)> an d h denotes the value of the longitudinal magnetic field. The 
function 4>(cp|m) generates a number of interesting correlation functions. We shall 
consider the two-point function 



(S\Sl) Tjh = \{2D 2 m d% -4D m dy+ l)*(q>K 
and the so-called emptiness formation probability 



(3) 

9=0 



P T (m) = lk4(9|ffl) = (II^Dr,* < 4 > 

7=1 

which is the probability to find a ferromagnetic string of length m in the equilibrium 
state with temperature T and magnetic field h, and is a special matrix element of the 
density matrix. D m in © denotes the 'lattice derivative' defined on any complex se- 
quence {z n )nen by D m z m =z m -z m -\, and e\ in © is the matrix 

We found in our previous works [9, 10] that the same auxiliary function a that 
determines the free energy of the model also contains all the information about the 
state of thermal equilibrium that is necessary to express the correlation functions as 
multiple integrals. As will become particularly clear in the XX limit, the function 
1/(1 + a) generalizes the Fermi function to the interacting case. Since we expect 
the Fermi function to appear naturally in all expressions for correlation functions, we 
think that the formulation of the thermodynamics of integrable models as developed 
in [22, 23] is particularly useful in the present context. 

The auxiliary function is the unique solution of the non-linear integral equation 

lna(3v) _ h 27sh 2 (T 1 ) r d(D sh(2rQln(l + a(G))) (5) 



T Tsh(X)sh(X + T\) ,/c 2ra sh(?i-co + ri)sh(X-a)-ri) 

Here the integration contour C is a rectangle in the complex plane. We have sketched 
it in figure[T] Like the inhomogeneity and the kernel in © it depends on the parameter 
T|, which is the only essential parameter of the model and parametrizes the anisotropy 
A in ([TJ as A = ch(r|). Knowing the auxiliary function on C we can calculate the free 
energy per lattice site according to the formula [9] 

h f d(0 sh(ri)ln(l + q((D)) 
n ' } 2 Jclni sh(co)sh(Q) + Ti) ' K ) 
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Figure 1 : The canonical contour C in the off-critical regime A > 1 (left panel) and in 
the critical regime |A| < 1 (right panel). For A > 1 the contour is a rectangle with sides 
at ±i| and ±y, where y is slightly smaller than ^. For |A| < 1 the contour surrounds 

the real axis at a distance |yj slightly smaller than ^ if A > and slightly smaller than 
|_ N if A<0 _ 



We can now state the main result of [9]. The generating function <J>(cp|m) of the 
correlation functions of the XXZ chain in the thermodynamic limit, as defined in ©, 
has the following integral representation: 



4>(<p|m) = £ 



m p (m— rc)(p 



iusf 



^,ysh(^-r|) 



to W LjLVr23riV sh (^-) 

' A sh(co y --^-ri) 
ill sh(Cy-^-Tl) 

where a(X) is the auxiliary function ©, where 



sh(co / ) 



2ra(l + o(co ; )) Vsh(o) ; '-ri) 



det„ Af (©;, C/t) det„ G(u) y , , (7) 



xmI r\ /r , A sh(ffl r ^-1l) " sh(tOy - ^ 

M ( »,, w = .(&,«.;) n uip^ +'(»;. ji sh( M ;- mf+ ^ 



f(£,GD) 



sh(ri) 



sh(^-co) sh(^-co + ri) ' 
and where the function G(X, Q is the solution of the linear integral equation 



(8) 
(9) 



sh(2Ti)G((0,C) 



27ii(l + a((o)) sh(A,-(o + r|)sh(A,-(fl-r|) 



(10) 



The contour C in the latter equation is the same as in the non-linear integral equation 
© (see figure[T]). T is any simple closed contour which lies inside C and encircles the 
origin counterclockwise. 
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An alternative expression for <3>((p|m) which involves the function a = l/o and 
which was also obtained in [9] is 



4>((p|m) = £ 



(-1)" 



n2 



fr [d^ / sh(Cy + Tl) yY ^ / sh(0) J ) 

JLi ir 2m V sh(Cy) ; ic 2ra( 1 + a((o ; ) ) \ sh(a>; + tj) 
sh((o y -^ + ri) 



det„ Af(a>,-, det ;7 G(co 7 , . (11) 



For the emptiness formation probability © we obtained several alternative expres- 
sions [8,9]. In the present context it will turn out to be convenient to work with the 
inhomogeneous, symmetric version 



1 



ml 



n 



J^Jc 2ra(l + a((o 7 ))_ 

Sh(0);-^) sh(O)y-^-Tl) 



sh((o 7 -03 fc -'ri) 



det m tfe k , (Oj) det m G((Oj£ k 

n^ =1 sh(^-^) 



(12) 



which follows from equation (1 1 1) of [9] in the limit (p — > —o°, and to perform the limit 
t,j — > 0, j = 1 , . . . , m, at a later stage. 

Formulae for the magnetization are the simplest realizations of the above integral 
representations, e.g. 



m(T,h) = (S^ T ,h = limP T (U) - \ = . 
J c,— »0 

Using ©, (fTTT) in the last expression on the right hand side we obtain [9] 



(13) 



m(T,h) 



1 

2 



Ju) G(co,0) 



1 



+ 



2m 1 + o(co) 2 y c 2jii 1 + a(to) 



dco G(u),0) 



(14) 



In the following sections we shall consider the limits r\ — > in/2 and r\ — > °° with 
JA — c fixed. In the former case A goes to zero and the Hamiltonian (UJ) becomes the 
XX Hamiltonian. This limit is also called the free Fermion limit. The latter limit is 
called the Ising limit, since in this limit the Hamiltonian © turns into the Hamiltonian 
of the one-dimensional Ising model. 



3 The XX limit 



In the XX limit r\ = irc/2, A = ch(r|) = we have 



sh(X±r|) =±ich(X), 



(15) 
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and the equations (l5l)- (fT2l simplify considerably. In particular, the kernel of the inte- 
gral equations © and (fTOb vanishes identically, and the auxiliary function a and the 
density function G are given by the corresponding inhomogeneities, 

»W = «p{4( A+ iipo)}' (16) 
G ^ = -M2^Q Y <17) 

Note that the function G(X, Q does not depend on the temperature and therefore agrees 
with its zero temperature limit. For this fact the following calculations bear many 
similarities with the zero temperature case studied in [17]. 



Magnetization 

We wish to transform all integrals over C into integrals over the 'Brillouin zone' 
[—71,71], To begin with a simple example, we consider the equation (TbH) for the mag- 
netization. Inserting ( fTBT ) and (fTTl) we obtain 

1 f d(0 1 

m(T,h) = ~- ; — t^^. (18) 

2 ic^ish(2( )l + exp(I(/ J + 4^)) 

The contour of integration C consists of a part C+ above the real axis and a part C- 
below the real axis (see figured])- Since the integrand in (fTSl) is analytic in the strips 
— 7t/2 < lm(0< and 0<ImG)<7t/2we may deform C± into contours parallel to the 
real axis passing through the points ±i7t/4. Then the parameterization 

1 in 

C+: [f,Jt]u[-7t,-f] ^C, (0( J p) = -arth(sin(p)) + -, (19a) 

C_:[-l,f]^C, a>(p) = -arth(sin(p))-- (19b) 

of C is uniform, and for every point raonC there is a unique p E [—71, jc] such that 

sh(2a>) = , ch(2co) = -itg (p) (20) 

icos(p) 

It follows the familiar expression 

m(T,h) = -- r ^ l —r, . (21) 

1 ' 2 y_ rt 27U + exp(^£2^) 



7 



Emptiness formation probability 



The formula for the emptiness formation probability is obtained in a similar way. We 
first insert dT31)-(IT7b into (fT2l and use the product formula for the Cauchy determinant. 
Then 



i 



m! 



n^ = iSh(03 ; --^)ch(0) ; --^) 



(22) 



and the homogeneous limit — > is obvious, 

m(m— 1) 



P T (m) 



m! 



jEJ^ri(l + a(a> y -))J n7=ish OT (2a) 7 



Ili<;<fc< m 4sh 2 ((0 J --co^ 



(23) 



Again the contour and the integrand are uniformly parameterized by (fl9t . in particular, 
using (|2H . 



1 



1 



1 



sin 



2 m 



sh (( 0j -( ,) = -ch(2a) J )ch(2( ,)--sh(2a) 7 )sh(2a),)-- co<p . )co<pk] 
for a// £ [ — 7C.7C] . Inserting (fT9l) . d20b and d24h into d23t we end up with 



(24) 



fr 

27! 1 + 



1 



4/ cos pi— h 

exp( j- 1 — 



J] 4sin 2 (^-^). (25) 

\<j<k<m 1 



This can be further simplified, if we take into account that rii<j<yt<m4sin ( — 
can be expressed as a product of two van-der-Monde determinants. Then 

* n dp e Kj-k)p 



P T (m) =det m 



-« 27t l+exp( 47c y- /i 



(26) 



j,k=l,...,m 



which is the known result, first obtained in [25]. For a careful discussion of the large 
distance asymptotics of P]{m) at zero and non-zero temperature see also [25]. 

Two-point function 

For the calculation of the two-point function (S\Sfy T , we use the integral representa- 
tion (fTTT) . Inserting (IT3T>-(IT7b and applying the product formula for the Cauchy deter- 
minant twice we find 



4>((p|m) = 
« ( e <P-l)" 



d(Hj f th(a)j 



^Vr 27ii 2ra(l + a(co_/)) \ tn (C 



det n 



sh(co J -^)_ 



(27) 
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where a = 1/a is given explicitly in (fT6b . The right hand side has a rather similar 
structure as in the zero temperature limit, and, at first instance, the same tricks apply 
(cf. [17]). Since the parameter (p enters (l27l) in the combination (e* — 1)'", only the 
three lowest order terms in the sum remain after applying the differential-difference 
operator at the right hand side of ©, and 



[S\S' m ) Th = \D 2 m {h+h)-D m h + \ 



4 ' 



where 
h 

h 



dt, 



doo 



/th(co) 



r2%iJ c 2Tii(l + a((D)) \th(Q J sh 2 (o)-Q 



n 



d&i 



f = \ Jr 2ni J c 2ra(l + a(tOy)) \ th(Cj 



th((B;) 



det? 



(28) 



(29a) 



(29b) 



The T-integrals in d29l can be calculated by means of residue calculus. For this 
purpose we deform the contour T as sketched in figure |2| Then 



dt, /th(a>)\ 



r 2%i 



i 



d% Ah(co)Y 



1 



sh 2 (o)-Q h2*x\HQ) sh 2 ((o-0 

1 



res 



VKQJ 



2m 



sh 2 ((O-0J^ C0 sh(2(D) 



(30) 



where the integral over f is zero, because the integrand is periodic with period in and 
vanishes asymptotically for Re C, — > ±°°. Inserting (QUI) into d29al ) and comparing with 
( fTSl we conclude that 



h=m[\-m{T,h)]=m[\-{^ T A 



(3D 



Similarly 
/ 2 



If- 



■ 2 

n 

■7=1 



II- 



2 

n 



J(Oj th m (tOy) 

2ra(l +o((Dy-)). 

J(0, 



,/r 2ra th m (Q sh(^ - ©i) sh(^ - ro 2 ) 



2m(l + a((Dy)). 



_sh((0i —0)2) 



th(coi) 
_th(a) 2 ). 



th(oo 2 ) 
_th(a>i)_ 



(32) 



Inserting the above expressions for I\ and 7 2 into (|2"81 and applying the difference 
operators we arrive at 



d(D d a th m ((0) 



27li(l + a((C0)) m 



a m cth ra ((o) 



27ii(l + a(co)) m 



(33) 
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Im 




" 2 



oo 



Re 



Figure 2: In order to calculate the integrals over C, in I\ and h the contour Y is deformed 
into f plus a small circle surrounding the pole at CO. 



This is already a rather simple expression involving only products of single integrals. 
We wrote it in a form that is particularly convenient for transforming it into an integral 
over the Brillouin zone: Since (see (1201 ) 



th(co) 



1 



ch(2co) - 1 



cth(co) sh(2oo) 



-le 



(34) 



we see at once that 



\ S l S m+l)T,h 



\ S l)T,h( S m+l)T,h - 



% dp cos(mp) 

Ott 1 i I h— 47cosp\ 
-7t 1% 1 +exp( j — L ) 



(35) 



for all m E N. 

Let us finally consider the zero temperature limit. If h is larger than the critical 
field 47, all spins are ferromagnetically aligned at zero temperature and the correlation 
function (l35l) vanishes. For \h\ <AJ the integral (l35l) in the zero temperature limit is re- 
stricted to the interval between two Fermi points ±arccos(^y), where the exponential 
function in the integrand vanishes, and 



( S l S m+l)T,h ( S l)T,h\ S m+l) 



T,h 



T 2m(jj) ~ 1 

2% 2 m 2 



(36) 



with Tj(x) being the jth Chebychev polynomial. Equations (031) and d3*6T ) are equivalent 
to formulae first obtained in [24] (for h = 0) and [1]. The asymptotics of the two-point 
functions was studied in [1]. 
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4 The Ising limit 

The XXZ Hamiltonian CQ) turns into the Hamiltonian of the one-dimensional Ising 
model if we set 7 = -| = c > 0, and take the limit T| — > °°. In this section we 
shall reproduce the known formula for the emptiness formation probability in the Ising 
limit [27] from the multiple integral representation (fT2l . As a prerequisite we need to 
calculate the auxiliary function a(X) and the inhomogeneous density function G(X, Q 
from their integral equations © and (fTOt . We recall that, in the off-critical regime, the 
integral equations determine the functions a(X) and G(X, Q inside the strip \lmX\ < 5, 
|ReA,| < (1 - 8)5 for any fixed 8 with < 8 < 1 . Using that 

limcth(X±Ti) =±1, (37) 

uniformly for all X with |ReA,| < (1 — 5)rj, and the elementary formulae 

sh(ri) 



sh(X) sh(X + ri) 
sh(2T]) 



sh(X-r|)sh(X + r|) 
we see that the integral equations ©, (flOb turn into 

h 2c , , , x f do 



cth(X)-cth(X + Ti), (38) 
cth(^-Ti)-cth(^+Ti), (39) 



lna(X) = -£-^(cth(A,)-l) + / — ln(l+a(co)) 
11 Jc rci 



(40) 



G(U ) = c*K-X)-i-/^^) (41) 

Jc 711 1 + a((Q) 

These equations are now valid for all X E C with |ImA,| < |. In particular, we may 
choose for X any point on the left or right edge of the contour C (see figure CJ) and 
then shift the left edge to — °° and the right edge to +°°. This way we obtain from both 
integral equations, d40l) and (HTI) . two algebraic equations that determine the limits 
a± = lim| Re ^|^ ±00 a(?i) and G± = lim| Re ^|^ ±00 G(?i, Q, respectively. The equations 
for the limits of the auxiliary functions read 

h—2c 1 + d+ 2c 

a±=Ke-T^T : k= (42) 

1 + a_ 

This is easily turned into a quadratic equation for K whose roots are 

K = e-r+^ (- S h(—) ± JshH-)+eA . (43) 




Shifting the left and right integration contours to =F°° the constant K appears in (l40l) . 
and we end up with 

a(X)=Ke-T-T cth M. (44) 
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The sign in d43t is fixed to plus by taking into account that a(X) in the high temperature 
limit is identically equal to 1 for arbitrary value of the anisotropy parameter A (see [9]). 
The function G(k, Q is obtained in a similar way. The limits G± satisfy 

G+ G- G+ G 

+ - , G =— r-^— + - (45) 



l + o + l + o_ l + a+ l + a_ 

which follows from (BIT) . We infer that G± are independent of C, and that G — G + = 2. 
Hence we can introduce a new function, say (a), such that 

G± = (g) T 1. (46) 

Inserting the latter definition into (H31) . using (|4"2"1) . d4"3t and solving for (a), we obtain 
the explicit expression 

(o)= , Sh( ^ =• (47) 

By shifting the integration contour to =p°o for fixed A, in (I4TT) and using (l45t . d46l we 
eventually arrive at 

G(U) = cth(C-A-) + (a). (48) 
Having determined the auxiliary function a(X) and the function G(k,Q we can 
now calculate the magnetization as a first simple application. Shifting the integration 
contour to =f°° in the first equation (TBI) and making use of d45t and (l46b we obtain at 
once 

m(T,h) = Q, (49) 

i.e. (o) = (o Z j)T,h which explains our notation. Not totally unexpected our formulae 
d4"9l and d4*7t are in agreement with the well-known result [2]. 



Emptiness formation probability 

Starting point for the calculation of the emptiness formation probability is again (fT2l . 
After (i) taking into account the simplifications due to (l37t -(l39l). (ii) using the formula 



det, 



axj + bsk 



Xi-S k 



(a + b) 



m—\ 



n 



\<j<k<m 



(Xj-X k )(s k -Sj) 



(50) 



aY[xk + b\[s k 

. k=\ k=l 

in order to simplify det ra G((0/,^)> and (iii) setting s k = e~ 2 ^ and changing the inte 
gration variables to xj = e~ 2(0 j, we arrive at 



2m! 



n 



dxj 



1 



;=l J? 2nix .il+Kexp 



h _ 2c ]+Xj 
T T l-Xi 



EI ( X j- X kj 
l j,k=l 
)+k 



(l + (a)m =l x k + (l-(a))UT=^k 
n™ k=l (xjSk) 



■ (51) 
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Here we denoted the transformed contour by 1 . It consists of two concentric circles 
about the origin. The left edge at, say, — b of the original contour of integration is 
mapped to a large circle of radius e 2b , while the right edge at, say, a is mapped to a 
small circle of radius e~ 2a (see figure 13). As the radius of the large circle is increased 
the auxiliary function on the contour tends to its limiting value a_ = Kq~^ 1 ~ 2c '' t . 
On the small circle it approaches a + = Ke~~( h+2c >> T if the radius goes to zero. Since 
the remaining factors in the integrand are rational in the xj with only simple poles at 
xj = Sk, k = 1 , . . . , m, and at xj = 0, we may calculate the integral by first fixing a 
large and a small circle, such that all poles of the rational part of the integrand, except 
the poles at zero, lie in the annulus between the circles, then replacing the auxiliary 
function by a on the large circle and a + on the small circle and finally using the 
residue theorem. After carrying out n integration this way the right hand side of (BTT) 
takes the form 



1 



2m! (1 + a-) n 



n 

7=1 



dxj 



1 



| I [Xj-X k ) 

.j,k=\ 

tfk 



(i+(^»nr=T^+(i-(^))n m fe=1 ** 

k^jl,-Jn 

n7 = 7iT% =1 (xj-sk) 

k^jlj-Jn 



-n(l-<o»- 



n&L i n"= J ( s j e x k ) n" = I sj e 



nm-n 2 
k=\ X k 



i+a+ n m k=l uu( s k- s k) rrprrr k=i (*,-^) 



(52) 



which is most easily verified by induction over n. 
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For n = m numerous trivial cancellations among the products occur, and we obtain 



*rH{$}) 



1 



;i+a_)' 



1 - (a) a, 



in 



S( 



This further simplifies due to the identity 



(53) 



En 



Sf 



k=i£=l S( Sk 

e^k 

which holds for any set of mutually distinct complex numbers Sk, k= 1 , 
ing J54t and <05]», J46i) into <|53]> we end up with 



l + <o> 



2(1 + 0- 



\m— 1 



(54) 



, m. Insert- 



(55) 



where (a) and a- are explicitly shown in equations d4*7b and (|42b . d43l . Quite remark- 
ably, the right hand side of (l55t is completely independent of the inhomogeneities {^}. 

In order to express the emptiness formation probability in more physical terms let 
us calculate the free energy in the Ising limit from the integral ©. Proceeding as above 
we find 



f(T,h) 



h 

-Tin 



rin(i + a_; 



ch (2r) 



+ Wsh 



2T 



Ac 

+ e t 



(56) 



Comparison of (l55t with (l49b and d56t then yields 

/ 1 \ n f(Jfi) i nh 

P T (n + l) = (i+m(r,/i))e^+2T 



(57) 



in accordance with [27] and with the general formula for the asymptotics of the empti- 
ness formation probability at finite temperature [5]. 



5 Conclusions 

Our aim in this article was to demonstrate the usefulness of finite temperature multiple 
integral representations for correlation functions, as obtained in [8-10], by showing 
that a number of known results for the correlation functions of the XX chain and of 
the one-dimensional Ising model can be reproduced by means of our formulae ©, 
dTTb and (fT2l . Usually the correlation functions of the XX and of the Ising chain are 
calculated from rather different, very special devices. The correlation functions of the 
XX chain are obtained by mapping the chain onto a model of free Fermions [24], while 
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the calculation of the correlation functions of the one-dimensional Ising model relies 
on the transfer matrix formalism with a spectral parameter independent 2x2 transfer 
matrix [2]. Here, for the first time, we have obtained correlation functions for the 
XX chain and for the Ising chain from a common device, namely the multiple integral 
representation, which interpolates for all values of A between the XX limit A = and 
the Ising limit A — > oo. 

We would like to see this work in a line with [28], where a high-order high- 
temperature expansion for correlation functions of the isotropic chain was obtained 
starting from the multiple integrals, and with [7], where, again from the multiple in- 
tegrals, some of the two- and three-site functions were computed numerically. We 
hope to report on further applications of the finite temperature integral formulae in the 
future. 
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